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Lattice determinations of quark mass have made significant progress 
in the last few years. I will review recent advances in calculations of charm 
and bottom mass, which are near to achieving percent-level precision and 
with fully controlled systematics. Precise knowledge of these parameters 
is of particular interest for precision Higgs studies at future accelerators. 
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1 Introduction 


Quark masses are fundamental parameters entering into the definition of the Stan¬ 
dard Model. Within the Standard Model picture, quark masses arise from Yukawa 
interactions with the Higgs held, and direct measurements of the Higgs couplings at 
the LHC are consistent with Standard Model predictions. High-precision studies at 
future accelerators such as the ILC will measure couplings at the per mil level [1]. 
In order to test the SM at this level, and to constrain and potentially discriminate 
between models of new physics detectable at this level, it is imperative to determine 
the quark masses to a corresponding level of precision. 

In recent years, considerable progress has been made in lattice calculations of 
quark masses, with groups now quoting charm and bottom mass values at around the 
percent or few-percent level. This is due to increasingly realistic simulations, and new 
techniques. State-of-the-art simulations include dynamical u, d, s, and frequently c 
quarks, with pion masses reaching their physical values, and typically at several lattice 
spacings. This increased realism translates into increasingly accurate results, and 
with fewer systematic errors. In order to reliably determine quark masses at the sub¬ 
percent level, it is important to have a variety of calculational techniques/strategies 
available, along with independent determinations from different groups. 

The outline of the rest of this article is as follows: Sec. 2 briefly discusses quark 
mass parameters in a general context, and how they are determined in lattice QCD 
simulations. In Sec. 3 I will discuss recent progress in the charm mass determinations, 
focusing on a promising method using current-current correlators. Sec. 4 will look 
at strategies and results for bottom mass determinations, and Sec. 5 discusses the 
important roled played by mass ratios. Sec. 6 presents a summary and discusses 
future prospects for these calculations. 

2 Quark mass and LQCD 

Quark masses are scheme and scale dependent quantities and can be viewed as in¬ 
put parameters that, along with as, specify QCD at the Lagrangian level. These 
parameters must ultimately be determined from experiment, but because quarks are 
conhned into hadrons the connection is necessarily indirect. In the absence of lattice 
simulations, one must focus on experimentally measureable observables which are 
1) sensitive to quark masses and 2) can be reliably computed in perturbation the¬ 
ory. One set of observables satisfying these criteria are derived from the the R-ratio. 
Much effort has gone into calculation of the relevant perturbation series, which are 
now known to N^LO [2, 3, 4]. As will be discussed in Sec. 3.1, one promising way 
to calculate rriq for heavy quarks via lattice simulations uses the same perturbative 
calculations, but substitutes experimental data with data from LQCD simulations. 
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Lattice QCD simulations are well suited for mass determinations, since the mass 
parameters are simulation inputs controlled by the “experimenter”. By changing the 
input masses, one can directly measure the resultant change in physical observables. 
In a standard LQCD simulation, one tunes the input masses in order to reproduce 
the masses of some low-lying hadrons - one for each quark in the theory. In this way 
one obtains (typically very precise) bare quark masses, but in the particular lattice 
regularization one happens to be using. In order to make contact with a continuum- 
regularized determination such as the MS scheme, one needs an additional calculation 
of the lattice to MS matching factor. This can be found using lattice perturbation 
theory or via non-perturbative renormalization (NPR) techniques. The ratios of bare 
quark masses in a given regularization are however immediately useful, as they are 
equal to renormalized mass ratios (up to lattice artifacts). 


3 Charm quark mass 

3.1 Current-current correlator method. 


The current-current correlator method uses moments of Euclidean-time twopoint 
functions, 

G{t) = a® x) J5(0, 0)). (1) 

X 

Here J 5 = and arrioh is the bare quark mass parameter in lattice units. In 

formalisms with sufficient chiral symmetry, the current J 5 is absolutely normalized. 
The correlator G{t) is UV finite, so that 

C(t)cont = C(f)latt + C(a^) (t 7^ 0). (2) 


The correlators G(f)iatt are the same ones used to compute pseudoscalar masses and 
decay constants, in which case it is the large-t exponential tail of the correlator that 
is of interest. For the mass calculation it is the small-t short distance behavior that 
is extracted via time-moments of G{t), defined as: 


T 

G^n,latt latt • 

t=0 


(3) 


The time-moments Gn have also been computed to N^LO in perturbation the¬ 
ory [2, 3, 4]. For u > 4, 

9n{<^MS^ 9 ) 


G 


n,pert 




n—4 


(4) 


Here is the MS quark mass at the scale p. The basic strategy to extract the 

quark mass is to compare Gn.cont, the continuum extrapolated Gn,iatt values, with the 
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perturbative expressions Gn,pert in Eq. (4) (evaluated at a scale fi ~ mh), and from 
these determine best-fit values for and For example, computing the 

continuum limit of G 4 ,iatt with physically tuned input charm masses moc, one can 
obtain ay^{mc), and then use this value in Gg to obtain rricijnc)- 

The HPQCD collaboration carried out an analysis in [5] using reduced moments, 
Rn, which are simply related to the time-moments as 


= Gi/Gf 


(5) 

( 6 ) 


Rn = — (G, 


(GJGW)1/("-") (n > 6). 


rriQc 


where Gn'^ are the tree-level results for the moments. Dividing by Gn'^ has the 


advantage of reducing lattice-spacing effects. In continuum perturbation theory. 


Ri — '^4(q^ms> ^^) 


(7) 

( 8 ) 



Here are the perturbative expressions given by appropriate powers of gnlgn \ with 
gn'^ the lowest order perturbative result. For a given one computes the values 


of Rn from Eq. (6) and gets an estimate of mdSmh) = Rn/rn{Smh), via Eq. (8) 


(the scale fi was taken to be 3mh). In this way the scale dependence of m^®(/r) is 


determined. 

The running of mdfi) was calculated this way in [5] using Uf = 2 + 1 + 1 HISQ 
ensembles. The n = 4, 6, 8,10 moments were computed using three different lattice 
spacings a ~ 0.12,0.09,0.06 fm and for seven input masses from rrih = rUc - O.Trufe. 


The extractions of mc(3m/i) from each of these data points are shown in Fig. 1 (left). 


along with the perturbative running. Fig. 1 (right) shows the corresponding estimate 


of as extracted from this data, ran to Mz and compared with results based on other 


experimental inputs. 

Estimates of rndfi) from time-moments are subject to a number of systematic 
uncertainties. The truncation of perturbation theory of course limits the precision. 
Fortunately the expansions of = 1 -|- are known for j = 1,2,3 and 

n < 10. The lattice moments are sensitive to condensate terms not captured in the 
perturbative expansions. These effects are suppressed like (AQCD/2m/i)'^, but they 
also grow with n. The lattice data also has cut-off effects, which grow like as {arrihY 
and decrease with increasing n; these trends are visible in Fig. 1. 


Fitting the moments data for n = 4, 6, 8,10 to Eqs. (7) and (8), HPQCD End 


mf®(3 GeV,n/ = 4) = 0.9851(63) GeV 
af^(3 GeV,n/ = 4) = 0.2545(37). 
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(9) 

( 10 ) 
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Figure 1: (Left) Data from [5] showing = Snih) extracted from lattice data 

and perturbation theory for moments n = 6 , 8 ,10 using Eq. ( 8 ). The green/blue/red 
data points correspond to lattice spacings of 0.12/0.09/0.06 fm. The gray band shows 
the evolution of the best-£t value for rric using perturbation theory. (Right) Value 
of a^^{Mz) from [5] compared with determinations based on various experimental 
inputs and a world average. 


These are compatible with earlier Uf = 2 + 1 results [ 6 ]. The compatibility of 
n + / = 2 + 1 and n/ = 2 + l + l results suggests that the effect of charm quarks in 
the sea can be treated perturbatively, to this level of precision. 

The JLQCD collaboration has recently utilized the current-current correlator 
method with nf = 2 + 1 domain-wall fermions to determine and [7]. Their cal¬ 
culation uses three lattices spacings a ~ 0.08, 0.055, 0.044 fm, and focuses on Rq, Rg, 
and Rio, from which they End 

mf®(3 GeV, Uf = 3) = 0.9936(91) GeV (11) 

af®(3 GeV, n/ = 3) = 0.2526(92). ( 12 ) 

3.2 Comparison of results 

In [ 8 ] the ETMG collaboration use lattice RI/MOM techniques to determine a mass 
renormalization factor 1/a) connecting the bare mass to the Rl-scheme mass, 

rrico, which is converted to the MS scheme using continuum 
perturbation theory. Unlike the current-current correlator method, which uses a 
heavy input mass to set the scale +, the RI/MOM calculation is extrapolated to 
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Figure 2: Comparison plot for determinations of = 4), computed from 

uj = 2 + 1 and n/ = 2 + l + l simulations. 


the chiral limit, and ETMC have generated mass degenerate rij = 4 ensembles for 
this purpose. The xQCD collaboration have also used RI/MOM methods for their 
uj = 2 + 1 determination [9, 10]. 

A comparison of recent lattice results for is shown in Fig. 2. 


4 Bottom mass 

It is challenging to directly simulate the b mass in relativistic lattice simulations, since 
one would like arribo 1 to keep discretization effects under control. Instead effective 
theories may be employed such as non-relativistic QCD (NRQCD) or heavy-quark 
effective theory (HQET). It has recently become possible with improved relativistic 
actions to approach the b mass, making extrapolation methods viable. 

4.1 NRQCD approach 

The NRQCD Hamiltonian is written as an expansion in v^, where u is a typical 
velocity of a 6 quark in the system of interest. For example, ~ 0.1 in the T meson. 
NRQCD calculations should be carried out with arribo > 1. This has the advantage 
that the b can be simulated using relatively coarse lattices, on the other hand it is less 
straightforward to extract continuum physics as compared to relativistic calculations. 
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The NRQCD current-correlator approach [11] is similar to the relativistic approach 
described in Sec. 3.1. One studies the time-moments of Euclidean-time two-point 
correlators. Unlike in the relativistic case, here the currents need to be normalized, 

JNRQCD _ jcont _ ^^3) 

Then the time moments are related to continuum perturbation theory, 

^NRQCD _ y2 drijctusy h) 

Constructing ratios of successive moments, the factors of Zy can be canceled. 
Because the continuum limit cannot be approached directly one instead studies rrih 
as a function of the moment number. Compared to the charm case, condensate 
contributions which grow with moment number are more suppressed at the heavier 
quark mass. A “plateau” in rub as a function of moment number implies that n is 
sufficiently large for discretization effects to be small. Such a plateau from [11] is 
shown in Fig. 3 (left). 

Results at three lattice spacings and with two different light-quark masses for 
n = 18 are shown in Fig. 3 (right). A £t to this data, including systematic errors, 
and perturbatively evolved to rub gives 

^/ = 5) = 4.196(23) GeV. (15) 

This result is compared with others in Fig. 4. It is signihcant that the values in 
the figure are calculated using a range of techniques. In [6] results are extrapolated to 
rrib from below, using a relativistic action as described in Sec. 3.1. This calculation is 
based on a different range of moment numbers, and uses a different action than [11]. 
The work of [12] uses the binding energy of T and Bg mesons, computed using NRQCD 
and lattice perturbation theory, to determine the heavy quark pole mass, which is 
then converted to the MS mass with continuum perturbation theory. 


(14) 


4.2 Ratio method 

The ETMC collaboration have used the ratio method [13] to extrapolate relativistic 
nj- = 2-|-l-|-l simulation results around the charm mass to the bottom mass [14]. 
The method is based on the expectation from HQET that 


m 


M, 


hi 


lim , 

pole 


= constant, 


( 16 ) 


where M^i is the mass of a heavy-light meson and is the heavy quark pole mass. 
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Figure 3: (Left) extracted from the moments of NRQCD current-current cor¬ 
relators at two different lattice spacings from [11]. (Right) Results from the n = 18 
moment as a function of lattice spacing and for two different light-quark masses. The 
gray band gives the continuum determination with the total error budget. 


They use simulation data consisting of ratios of meson masses, Mhi{mh) / Mhiiruh/ A), 
computed for a series of masses rrih around the charm mass, e.g.: = rUc, 

= Xrric, ..., = X^rric- These ratios have the advantage that discretization 

effects proportional to {arrihY cire largely canceled. From this data they construct 
the function 


y{mh,X) 


Mhi{mh) p{mh/X) 

Mhi{mh/X) p{mh) 


(17) 


The functions p{mh) on the r.h.s. of Eq. (17) relate the pole mass to the MS mass and 
are known to N^LO in perturbation theory. y{mh, X) satishes limm^->.oo y{^h, A) = 1 
on account of Eq. (16), and so its value can interpolated between the charm region 
and the static limit using a motivated ht ansatz. Rewriting Eq. (17), the combina¬ 
tion Xy{mh, fhen a known function that evolves Mhi{mh/X) to Mhi{mh). 

Choosing A such that Mhi{rri/^) = for some N, they determine the b mass 

from = X^rric- 


5 Mass ratios 

Mass parameters are inputs to lattice QCD simulations, these are pure numbers 
(amo) corresponding to masses expressed in units of the lattice spacing a. There 
is one bare mass parameter for each quark in the simulation, and these must be 
tuned to reproduce the physics of QCD. The bare mass parameters are tuned by 
measuring low-energy observables such as meson masses, and requiring that these are 
be equal to their physical values. After this set of observables has been used to tune 
the simulation parameters, one has a set of numbers {{amudo)i {arriso), {arrico)}. The 
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Figure 4: Comparison plot for determinations of = 5), computed from 

uj = 2 + 1 and n/ = 2 + l + l simulations. 


bare lattice inputs are defined at the cutoff scale and depend on the details of the 
discretization. However, ratios of input masses are equal to the ratios of MS masses, 
up to discretization effects that vanish in the continuum, 

amlp ^ ^ ^ 

am2o 



Thus once the MS mass is known for one quark in the theory, this can be converted 
to the MS masses for the others using the input mass parameters. 

An example of this is shown in Fig. 5 (left), for the input ratio mpc/^os from [5]. 
In the continuum HPQCD hnd that 

= 11.652(65) (19) 

Using their result for m^®(/i) from the current-current correlator method discussed 
in Sec. 3.1, they obtain 

mfS (3 GeV,n/ = 3) = 84.7(7) MeV, (20) 

which is the most precise estimate to date. Fig. 5 (right) shows a result from [5] using 
input mass ratios to obtain Here the input mass is increased from ruco towards 









m^o, and finally an extrapolation performed to obtain 


mb(/J,,nf) 

{rrih, nf = 5) 



{anioc)^ 


4.528(54) (21) 

: 4.162(48) GeV (22) 



Fignre 5: (Left) Continnnm extrapolation of the bare qnark mass ratio from [5]. 
(Right) Extrapolation of in to nsing simnlation data from [5]. Magen- 
ta/green/blne/red points are from lattice spacings of 0.15/0.12/0.09/0.06 fm. 


6 Conclusion 

Recent progress in lattice determinations of charm and bottom qnark mass was re¬ 
viewed. In order to achieve (snb-)percent level nncertainties for these qnantities, it is 
important that determinations come both from a variety of calcnlational strategies, 
and via independent measnrements from different gronps. 

The most precise qnoted valnes for c mass presently come from calculations of 
current-current correlators, comparing these to perturbation theory, where a heavy 
(~ rric) input mass sets the scale fi. The precision in the value of the charm mass can 
be cascaded to the other masses using bare quark-mass ratios, which are determined 
in the tuning of simulation parameters to their physical values. 

Calculations of b mass are done either using an effective-theory framework for the 
b quark or extrapolating relativistic simulations from lower-mass region where dis¬ 
cretization effects are under control. Extrapolation methods will continue to improve 
as ensembles with smaller lattice spacings become available. First steps have been 
taken towards a fully relativistic treatment of the b quark [5, 6]. This will lead not 
only to more precise values for the b itself, but through the use of mass ratios should 
improve determinations of the other quark masses as well. 
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